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ABSTRACT 


Multiple  coherence  gives  a  quantitative  measure  versus 
frequency  of  how  well  a  linear  combination  of  n  input  channels 
can  match  the  (n  +  1) st  channel  in  a  seismic  array.  If  the 
inputs  can  match  the  output  exactly,  then  the  multiple  coherence 
is  unity  and  only  n  channels  for  short  period  noise  fields  at 
UBSO  and  TFSO„ 

The  multiple  coherence  of  the  noise  at  UBSO  and  TFSO  short 
period,  vertical  component  arrays  is  high  (greater  than  0.9)  over 
the  microseismic  frequency  band. 

The  decay  of  the  multiple  coherence  of  the  noise  with  in¬ 
creasing  frequency  is  faster  at  TFSO  than  at  UBSO  and  faster 
at  UBSO  than  at  LASA. 


1.  INTRODUCTION 

Most  basic  data  processing  techniques  for  signal  enhancement 
or  identification  depend  upon  the  structure  of  the  noise  within 
the  seismic  array.  If  some  of  the  coherent  noise  is  due  to  site 
characteristics  such  as  consistently  coherent  noises  from  par¬ 
ticular  direction,  then  techniques  using  multiple  coherence  will 
help  to  isolate  these  consistent  linear  relations.  Many  optimum 
filters  for  estimating  the  signal  take  account  of  these  linear 
relations  implicitly  by  weighting  with  the  inverse  of  the  spectral 
noise  matrix.  However,  one  cannot  tell  whether  the  coherent  noise 
involved  is  due  to  noise  generating  events  which  cannot  be  pre¬ 
dicted  or  controlled.  Thus  the  filters  must  be  recalculated  over 
a  period  of  noise  recording  immediately  prior  to  the  arrival  of 
each  single  signal.  Part  of  the  coherent  noise  generated  within 
the  array  may  be  due  to  various  causal  factors  for  a  particular 
array.  If  so,  we  can  learn  something  about  these  factors  by 
examining  the  linear  relations  between  the  various  array  elements. 
A  potential  benefit  here  is  that  a  consistent  linear  model  rela¬ 
ting  the  different  sub-elements  would  eliminate  the  need  for  com¬ 
puting  a  different  set  of  filter  coefficients  for  each  event. 

The  multiple  coherence  function  can  indicate  how  many  seis¬ 
mometer  outputs  in  an  array  are  necessary  to  properly  determine 
the  seismic  noise  field.  If  there  are  n  independent  seismic  noise 
components,  then  the  multiple  coherence  function  would  be  unity 
when  (n  +  1) st  seismometers  are  placed  in  an  array  to  measure 
seismic  noise  records.  If  part  of  the  background  is  composed  of 
incoherent  noise,  then  the  multiple  coherence  function  would  in¬ 
dicate  the  percentage  of  coherent  noise  present  and  the  number 
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of  seismometers  necessary  to  define  this  coherent  noise.  The 
filter  relations  determined  by  the  multiple  coherence  computa¬ 
tions  can  then  be  used  in  array  summation  to  bring  the  noise  into 
destructive  interference*. 

This  analysis  does  not  guarantee  that  such  optimum  proces¬ 
sing  is  possible.  For  example,  if  the  noise  and  signal  propaga¬ 
tion  characteristics  across  the  array  are  identical,  no  velocity 
filtering  scheme  can  be  expected  to  separate  the  two  ever  though 
the  multiple  coherence  might  be  unity. 

The  multiple  coherence  function  is  the  frequency  domain 
equivalent  of  the  prediction  error  filter  in  time.  If  n  input 
seismic  traces  predict  the  (n  +  1) st  trace  in  an  array  completely, 
then  the  multiple  coherence  will  be  unity  and  a  prediction  error 
filter  could  be  used  to  exactly  predict  this  (n  +  1) st  output. 

In  fact,  linear  filter  relations  derived  by  the  multiple  coherence 
program  produce  an  estimate  of  the  (n  +  1) st  trace  which,  when 
subtracted  from  the  actual  (n  +  1) st  trace,  given  a  prediction 
error  trace.  Thus  the  combination  of  the  filter  derived  in  the 
multiple  coherence  program  and  the  subtraction  operation  produces 
a  prediction  error  filter  as  shown  in  the  following  diagram. 


*For  the  mathematical  description  of  the  multiple  coherence  com¬ 
putation,  See  Appendix  1. 
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Prediction  Error  Filter 


The  first  objective  of  this  study  is  to  use  the  multiple 
coherence  function  to  estimate  the  degree  of  predictability  of 
the  short  period  noise  field  at  UBSO  and  TFSCb  The  multiple 
coherence  from  these  two  vertical  component  arrays  we  will  com¬ 
pare  with  the  similar  analysis  for  the  short  period  LASA  arrays 
(See  SDL  Report  #190,  "Multiple  Coherence  of  Short  Period  Noise 
at  LASA",  W.C.  Dean,  E,F„  Chiburis  26  June  1967)., 

The  second  objective  is  to  determine  from  multiple  coherences 
the  number  of  independent  components  comprising  a  given  noise  field 
and  the  percentage  of  incoherent  noise  which  cannot  be  cancelled  by 
any  kind  of  multichannel  filtering.. 

The  third  objective  is  to  determine  the  stationarity  proper¬ 
ties  of  the  noise  fields  We  accomplish  this  by  applying  the  mul¬ 
tiple  coherence  program  to  three  different  time  samples  from  the 
same  arrays  Then  the  multiple  coherency  filters  derived  from  the 
first  time  sample  are  applied  to  the  other  two  time  samples.  If 
the  filters  derived  in  the  first  time  sample  have  done  a  good  job 
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of  predicting  the  noise  field  in  all  three  samples,  then  the 
data  are  said  to  be  stationary.  On  the  other  hand,  if  the 
filters  from  the  first  time  sample  do  a  progressively  poorer 
job  of  predicting  the  noise  in  the  other  samples  relative  to 
the  filters  associated  with  those  samples,  then  the  noise  is 
non-stationary  to  some  degree.  This  deterioration  in  predict¬ 
ability  of  the  multiple  coherence  filter  can  quantitatively 
measure  the  non-stationarity  of  the  data**. 


**  For  a  theoretical  discussion  of  the  stationarity  computation. 
See  Appendix  II. 


4- 


2.  DESCRIPTION  OF  DATA 

We  computed  the  multiple  coherences  from  arrays  of  short 
period,  vertical  component  seismometers  at  UBSO  and  TFSO.  The 
output  seismometer  was  the  center  element,  seismometer  Z-10  at 
UBSO  and  seismometer  Z-21  at  TFSO.  The  array  diagrams  appear 
on  the  figures  showing  the  multiple  coherence  plots.  At  each 
observatory  the  first  inputs  added  were  from  the  outermost 
elements  of  the  array.  The  number  of  input  channels  increases 
from  2  to  a  maximum  of  9.  The  multiple  coherence  is  plotted  vs. 
frequency  with  the  number  of  input  channels  as  the  parameter. 
Since  adding  one  additional  channel  to  the  set  of  inputs  may 
not  increase  but  can  never  decrease  the  amount  of  input  infor¬ 
mation  about  a  noise  field,  the  multiple  coherence  must  be  a 
monotonically  increasing  function  with  increasing  number  of 
inputs. 

The  samples  were  slightly  over  three  minutes  long.  The 
noise  examined  at  each  array  came  from  two  samples,  one  recorded 
on  August  1966  and  the  second  in  January  1967.  Stationarity 
tests  were  conducted  on  the  August  samples  where  three  succes¬ 
sive  3-minute  samples  were  recorded. 
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3.  RESULTS 


Multiple  Coherence 

Figure  1  shows  the  multiple  coherences  versus  frequency 
for  two  time  samples  separated  by  three  minutes  and  recorded  at 
UBSO  on  25  August  1966.  A  diagram  of  the  array  elements  chosen 
is  shown  a^  the  center  of  the  figure.  The  ordering  of  the  inputs 
is  from  the "outef most  to  the  closest  and  is  listed  in  the  figure. 
The  multiple  coherences  are  computed  every  .1  cps  over  a  range 
from  .1  to  2.6  cps.  The  number  of  points  in  each  sample  is  2000. 
The  number  of  lags  computed  in  the  correlation  function  is  50. 

The  data  sampling  rate  is  10  points  per  second. 

A  similar  plot  of  multiple  coherence  versus  frequency  for 
a  sample  recorded  on  7  January  1967  at  UBSO  is  shown  in  Figure  2. 
The  number  of  lags  used  in  the  correlation  is  sufficient  to  detect 
any  propagating  noise  component  with  a  velocity  greater  than  .5  km 
per  second.  Thus  lack  of  sufficient  lags  cannot  be  a  cause  of  low 
multiple  coherence.  The  multiple  coherences  are  high  for  the  fre¬ 
quencies  in  the  microseismic  range  but  fall  to  low  levels  for  fre¬ 
quencies  froia  1  to  2  cps.  As  more  inputs  are  added  the  multiple 
coherence  remains  high  for  higher  frequencies.  With  9  inputs,  the 
multiple  coherence  stays  above  .8  or  .9  for  all  frequencies  less 
than  1  cps.  However,  with  this  many  input  channels  the  closest 
inputs  are  only  H  km  from  the  output. 

Figures  3  and  4  show  similar  effects  at  TFSO.  Here  the 
multiple  coherence  falls  to  a  low  value  at  frequencies  well  below 
those  at  UBSO  However,  one  must  again  note  that  the  dimensions 
at  the  TFSO  array  are  greater  than  those  at  UBSO.  When  we  compare 
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the  multiple  coherences  at  the  two  arrays  with  the  minimum  seis¬ 
mometer  separations  from  the  output  trace  at  UBSO  held  about 
equal  to  those  at  TFSO,  we  still  find  that  the  multiple  coherence 
at  TFSO  decays  faster  with  increasing  frequency  than  that  at  UBSO 

In  comparing  the  multiple  coherences  at  UBSO  and  TFSO 
with  those  within  a  subarray  at  LASA  find  that  the  multiple 
coherences  at  all  three  arrays  are  large  for  the  microseismic 
frequencies  and  decay  with  increasing  frequency.  The  multiple 
coherences  at  LASA  remain  high,  greater  than  .8  or  .9,  out  to 
frequencies  higher  than  either  UBSO  and  TFSO.  The  closest  seis¬ 
mometers  at  LASA  have  a  spacing  similar  to  those  found  at  UBSO. 
Thus  wc  see  that  the  noise  field  at  LASA  is  more  coherent  with 
8  or  9  inputs  than  that  at  UBSO  or  TFSO  in  the  signal  band. 

Power  Spectra 

Spectra  for  two  of  the  multiple  coherence  examples  on 
Figures  1-4  are  shown  on  Figures  5  and  6.  The  figures  show  the 
output  power  spectra  and  the  range  of  input  spectra  for  8  inputs 
on  Figure  5  and  for  5  inputs  on  Figure  6.  These  two  spectra  are 
representative  of  all  spectra  for  the  multiple  coherence  examples 
computed  in  this  report. 

Stationarity  Tests 

The  multiple  coherence  program  derives  a  set  of  n  filters 
for  the  n  input  seismograms  which  together  provide  the  best  linear 
estimate  for  the  (n  +  1)  st  seismic  trace.  The  difference  between 
the  observed  (n  +  1) st  trace  and  the  best  estimate  is  the  error 
trace.  If  the  multiple  coherence  is  unity,  the  prediction  is 
perfect  and  the  error  trace  will  be  zero.  If  we  form  the  ratio 


of  the  error  spectra  over  the  observed  spectra,  we  can  get  a 
measure  of  the  reduction  in  noise  power  possible  from  the  theoret¬ 
ical  optimum  filters.  Thus  the  db  improvement  as  a  function  of 
frequency  can  be  expressed  as 

db  =  10  log  (error/observed)  . 

The  prediction  error  filters  will  do  the  best  job  in 
eliminating  the  noise  background  when  they  are  applied  to  the 
noise  sample  from  which  they  are  derived.  However,  if  the  noise  is 
stationary,  the  same  filter  could  be  expected  to  do  nearly  as  well 
when  applied  to  later  time  samples  from  the  same  array.  Figure  7 
shows  the  expected  noise  reduction  in  db  when  the  prediction  error 
filters  that  were  derived  from  the  first  time  sample  are  applied 
to  the  first  time  sample.  In  addition  these  same  filters  are 
applied  to  the  second  and  third  time  samples.  The  three  3-minute 
time  samples  are  adjacent  to  each  other.  Figure  7  shows  that  re¬ 
duction  of  the  noise  background  from  a  prediction  error  filter  is 
expected  izo  be  over  30  db  in  the  fitting  interval  for  the  micro- 
seismic  frequencies  and  within  2-3  db  of  this  level  on  adjacent 
time  samples.  At  higher  frequencies  the  expected  noise  reduction 
in  the  fitting  interval  falls  steadily  until  only  1  to  2  db  re¬ 
duction  is  expected  for  frequencies  higher  than  1.5  cps.  In  the 
adjacent  time  samples  the  expected  noise  reduction  is  within  1-2  db 
less  over  all  frequencies. 

Figures  8,9,  =*nd  10  show  similar  plots  over  expected  noise 
reduction  for  3  adjacent  time  samples  at  TF.30.  On  Figure  8,  sample 
1  is  the  fitting  interval.  On  Figure  9,  sample  2;  and  on  Figure  10 
sample  3  is  the  fitting  interval.  The  curves  at  TFSO  show  results 
similar  to  that  at  UBSO.  The  main  differences  are  that  the  noise 
reduction  in  the  microseismic  band  at  TFSO  is  a  bit  lower  than 
that  expected  at  UBSO  and  the  expected  noise  reduction  falls  to 
essentially  zero  db  at  1.2  cps  instead  of  1.5  cps  at  UBSO 
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4.  CONCLUSIONS 

1.  The  multiple  coherence  of  the  noise  at  UBSO  and  TPSO 
short  period,  vertical  component  arrays  is  high  (greater  than  0..9) 
over  the  microseismic  frequency. 

2.  The  decay  of  the  multiple  coherence  of  the  noise  with 
increasing  frequency  is  faster  at  TFSO  than  at  UBSO  and  faster 
at  UBSO  than  at  LASA  subarrays. 

3.  The  noise  has  a  low  multiple  coherence  (less  than  .4) 
with  up  to  9  inputs  at  UBSO  and  TFSO  for  frequencies  greater  than 
1.5  cps.  Thus  the  noise  is  largely  incoherent  over  the  upper  part 
of  the  signal  band  at  UBSO  and  TFSO. 

4.  The  noise  at  UBSO  and  TFSO  over  adjacent  time  samples  is 
stationary  only  over  the  microseismic  frequency  band. 
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FIGURES 


1.  Multiple  Coherences  vs.  Frequency  with  1-7  Inputs  for  tvo 
noise  samples  recorded  at  UBSO,  August  1966. 

2.  Multiple  Coherence  vs.  Frequency  with  2-8  Inputs  for  two 
noise  samples  recorded  at  UBSO,  January  1967. 

3.  Multiple  Coherences  vs.  Frequency  with  2-5  Inputs  for  two 
noise  samples  recorded  at  TFSO,  August  1966. 

4.  Multiple  Coherences  vs.  Frequency  with  2-9  Inputs  for  cwo 
noise  samples  recorded  at  TFSO,  January  1967. 

5.  Spectra  of  short  period  vertical  compoi.  its  recorded  at 
UBSO,  August  1966. 

6.  Spectra  of  short  period  vertical  components  recorded  at 
TFSO,  August  1966. 

7.  Expected  noise  reduction  from  a  prediction  error  filter 
applied  to  the  fitting  interval  and  two  adjacent  time 
intervals  for  noise  recorded  at  UBSO  in  August,  1966. 

8.  Expected  noise  reduction  from  a  prediction  error  filter 
applied  to  the  fitting  interval  and  two  adjacent  time 
intervals  for  noise  recorded  at  TFSO  in  August  1966. 

The  fitting  interval  is  the  first  interval. 

9.  Expected  noise  reduction  from  a  prediction  error  filter 
applied  to  the  fitting  interval  and  two  adjacent  time 
intervals  for  noise  recorded  at  TFSO  in  August  1966. 

The  fitting  interval  is  the  second  interval. 

Expected  noise  reduction  from  a  prediction  error  filter 
applied  to  the  fitting  interval  and  two  adjacent  time 
intervals  for  the  noise  recorded  at  TFSO  in  August  1966. 
The  fitting  interval  is  the  third  interval. 
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igure  1.  Multiple  Coherences  vs.  Frequency  with  1-7  Inputs 

for  two  noise  samples  recorded  at  UBSO,  August  1966 


Multiple  Coherence 


Figure  2  Multiple  Coherence  vs.  Frequency  with  2-8  Inputs 

for  two  noise  samples  recorded  at  UBSO,  January  1967. 
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Figure  3.  Multiple  Coherences  vs.  Frequency  with  2-5  Inputs 

for  two  noise  samples  recorded  at  TFSO,  August  1966. 
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Figure  5. 
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Figure  7. 


Expected  noise  reduction  from  a  prediction  error 
filter  applied  to  the  fitting  interval  and  two 

,arf2«C®nt  time  intervals  noise  recorder  at 
UBS0  m  August,  1966. 
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Figure  8.  Expected  noise  reduction  from  a  prediction  error 
filter  applied  to  the  fitting  interval  and  two 
adjacent  time  intervals  for  noise  recorded  at  TFSO 
in  August  19t6.  The  fitting  interval  is  the  first 
interval . 
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Figure  9.  Expected  noise  reduction  from  a  prediction  error 
filter  applied  to  the  fitting  interval  and  two 
adjacent  time  intervals  for  noise  recorded  at 
TFSO  in  August  1966.  The  fitting  interval  is  the 
second  interval. 


Figure  10.  Expected  noise  reduction  from  a  prediction  error 
filter  applied  to  the  fitting  interval  and  two 
adjacent  time  intervals  for  the  noise  recorded  at 
TFSO  in  August  1966.  The  fitting  interval  is  the 
third  interval. 


APPENDIX  I 

♦Multiple  Coherence  Functions 


Consider  a  collection  of  q  clearly  defined  inputs  x^(t); 
i  ■  1*2,...  ,q,  and  one  output  y(t),  as  pictured  in  Figure  5.12. 
Let  G£  (f)  =  Gii(f)  be  the 


Figure  5.12  Multiple-input  linear  system. 


power  spectral  density  function  for  x^t),  and  G^tf)  be  the 

cross-spectral  density  function  between  x. (t)  and  x.(t).  De- 
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fine  the  N  x  N  spectral  matrix  by 
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♦This  explanation  of  multiple  coherence  functions  was  taken  from 
"Measurement  and  Analysis  of  Random  Data",  Bendat,  J.  S.#  and 
Piersol,  A.  G.,  John  Wiley  and  Sons,  1966.  For  more  detailed  theo¬ 
retical  developments  and  discussions  of  multiple,  partial  and  mar¬ 
ginal  coherence  functions,  see  this  text. 


The  ordinary  coherence  function  between  x^(f)  and  x^(t)  is 


defined  by 


"  U)  ajU)  cm 


(2) 


The  multiple  coherence  function  between  x^(t)  and  all  other 
inputs  x^(t),  x^t)  ,...  .excluding  x^{t),  is  defined  by 

y, ,.*(/)  =  i  -  [CM)  CM)]'1  (3, 

where  G1(g)  denotes  the  ith  diagonal  element  of  the  inverse 
matrix  Gxx(f)  ^"associated  with  Eq.  (1)  .  The  ordinary  and 
multiple  coherence  functions  are  both  real-valued  quantities 
which  are  bounded  by  zero  and  unity.  That  is. 


o  <;  y,%f) 

o  <;  Yi'jj) 

The  multiple  coherence  function  is  ?  measure  of  the  linear 
relationship  between  the  time  history  at  one  point,  and  the  time 
histories  at  the  collection  of  other  points.  That  is,  the  mul¬ 
tiple  coherence  function  indicates  whether  or  not  the  data  at 
all  of  the  other  points  linearly  produce  the  results  at  a  given 
point . 


^  I 
£  1 


(4) 
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APPENDIX  2 


Theoretical  Development  of  The  Stationaritv  Relation a 

A.  Noise  Reduction  Within  The  Fitting  interval 

A  number  of  useful  statistical  measures  such  as  ordinary 
and  multiple  coherence  can  be  used  as  tools  to  indicate  the  amount 
of  noise  reduction  feasible  in  a  multiply  coherent  array.  The 
basic  linear  model  Which  determines  the  db  reduction  possible  in 
the  noise  field  by  multiple  coherence  filtering  relates  a  refer¬ 
ence  element  (trace)  y(t)  of  an  array  to  the  other  elements,  say 
x^(t),  Xjtt),  ....  Xp(t)  in  the  array  through  the  linear  model 

P  , 

y(t)  "  Z  J  xlt  <*-«>  dtt  (1) 

k-1  -• 

Generally  we  determine  h^t)  as  the  time  invariant  linear  filter 
that  makes  the  mean  square  error  between  y(t)  si!  its  predicted 
value  a  minimum,  i.  e. 

P  . 

E  I  y(t)  -  £  J  h^  (a)  x^ft-a)  da|2  -  min  (2) 

k-1  -• 

«Aiich,  by  the  usual  orthogonality  principle  (see  Papoulis  1), 
yields  the  condition 

P 

Ey(t)  x^(t+T)  “  £  J  E  x^  (t-a)  x^  (t+T )  d(t 
k-1 

t.  -  1,2,  p, <T<—  (3) 
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or  p 

Vt(T)  "I  J  Va)  \ x  <T+a>  <*a  (4) 

k-1  *  1 

which  by  taking  Fourier  transform*  implies  that 
P 
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k-1  *  4 

Now,  the  mean  square  error  can  be  written 

P  p 
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’*"*  a2<«)  i.  the  multiple  coherence  end  (1  -  o2(u,)  )  m.„ur.. 
the  reduction  in  power  poeeible  at  the  frequency  «.  with  e2(e)  -  1 

the  Man  equate  error  ie  aero  and  with  «2(e)  -  0  the  wean  equate 
error  is  just 

m 
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which  is  the  original  power  in  the  process  y(t).  Now  the  quantity 
(1-a  (id)  )  Syy(<»)  represents  the  residual  power  at  each  frequency 
after  the  best  linear  estimate  of  the  form  (1)  has  been  subtracted 
out.  Hence  the  db  reduction  in  power  at  each  frequency  is  just 
the  ratio  of  the  output  power  of  the  residual  (see  equation  (2)  ) 
to  the  input  power  in  y(t)  or 

S  «.) 

IH(u>)  -  10  log  -g  — j —  -  10  log  (1  -  a2  (id)  )  (8) 

yy 

2 

where  n  (u>)  is  the  multiple  coherence  and  S  (ou)  is  the  power 
spectrum  of  the  error  process 
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B-  Noise  Reduction  Outside  The  Fitting  Interval 

We  would  also  like  to  determine  the  noise  reduction  in  db 
which  would  result  from  using  a  set  of  filters  9k(t),  k-1,  ...,p 
which  have  been  derived  either  from  another  fitting  interval  pr 
from  theoretical  considerations.  To  accomplish  this  let  h]c(t) 
be  the  optimal  filters  for  the  time  under  investigation  and  let 
gk^  an^  ot^er  set  of  filters  whose  mean  square  error  is  to 
be  compared  with  h^t).  The  mean  square  error  of  the  g  filters 
can  be  written  using  the  orthogonality  principle  as 
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Hence,  if  we  call  the  new  error  /’(t)  we  have 

(t)  =  y(t)  -  £  j  gk(„)  Xk(t.„)  da 

k=l 


with  power  spectrum  S  (un  <-ha 

the  new  value  for 

in  the  Sk(t)  filters  tvould  be 
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Equation  (12)  shown  that  the 


expressed  in  terms  of  the 


improvement  in  the  g.  <t>  filters  i 


correction  term  which 


improvement  in  the  hk(t>  filters 


is  zero  when  H  =  g. 
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The  improvement  values  I„(ui)  and  I„(ui)  in  equations  (8)  and  (12) 

H  G 

are  those  shown  in  the  main  body  of  the  report. 
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